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^ ; Abstract 

^ , We present single-particle and thermodynamic properties of the half-filled 

single-band Hubbard model in 2D calculated in the self-consistent fluctuation 
exchange approximation. The low-energy excitations at moderate tempera- 
tures and small U are quasiparticles with a short lifetime. As the temperature 
is lowered, coupling to evolving spin fluctuations leads to the extinction of 
these quasiparticles, signaled by a weak pseudogap in the density of states 
and by a positive slope in Re S(ki;',e) and a local maximum in [Im S(ki?,e)| 

at e = 0. We explain these results using a simple spin-fluctuation model. 
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The 2D single-band Hubbard model plays a central role in efforts to understand the 
behavior of electrons near the Fermi surface in the cuprate superconductors and their parent 
compounds [|I|. The model is characterized by a nearest-neighbor hopping energy, t, and 
an on-site Coulomb energy, U. For a half-filled band, the T = state is believed to be an 
antiferromagnetic insulator for all f/ > [0,0. For T > 0, the Mermin- Wagner theorem 
precludes the existence of long-range AFM order in 2D, but with strong coupling {U ^ 8t) 
the low-temperature electronic state is almost certainly a Mott insulator. With U = 4t and 
O.lOt ^ T ^ 0.25t, conflicting results have been obtained from quantum Monte Carlo (QMC) 
calculations of the one-electron spectral weight function, y4(k, e), depending on lattice size 
and especially on the method used to extract A{\<.,e) from the Green's function G(k, r) 
produced directly by QMC. Spectral functions on the Fermi surface (FS) produced by the 
maximum entropy technique show a single peak whenever the lattice is larger than the 
AFM correlation length, but develop a pseudogap on smaller lattices [^]. In contrast, recent 
calculations using the method of singular value decomposition yield a pseudogap in both 
the spectral function on the FS and in the total density of states, N{e) = N^^ X]k^(k, 
even for lattices larger than the correlation length 

We report calculations of yl(k, e) and N{e) at half-filling and moderate U (< 4.8t) using 
the fluctuation exchange approximation (FEA), a self-consistent conserving approximation 
that has been applied to the 2D Hubbard model in a number of recent papers In 
particular, the FEA has been used to argue for a d-wave superconducting transition in the 
high-Tc cuprates ; the need to evaluate these claims adds to the importance of knowing 
what the FEA predicts (rightly or wrongly) for the normal state at half-filling. Compared 
to QMC, the FEA has the disadvantage that it is inherently an approximation, though 
imaginary-time Green's functions from the FEA and QMC agree surprisingly well at half- 
filling and moderate U 0. For studying the spectral function and DOS, the FEA has 
several important advantages over QMC: (1) There is no inherently statistical error in the 
FEA results, which removes most of the uncertainty in extracting A(k, e) from G'(k, r). (2) 
FEA calculations are possible for large enough lattices (typically 128 x 128) to ensure that 
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the profile of A(k, e) represents tlie infinite lattice limit, and does not in part reflect the small 
number of decay channels available to a low-energy quasiparticle on a small lattice. (3) The 
FEA can cover a wider range of temperatures than are currently accessible to QMC (4) 
The FEA provides real-frequency self-energies and fluctuation propagators that explain the 
origin of structures in A(k,e) and directly test the applicability of the Fermi liquid picture. 

We find that in the FEA a paramagnetic non-Fermi liquid state evolves with decreasing 
temperature. The spectral functions on the FS show no pseudogaps over the range of 
temperatures covered by QMC calculations, but are exceptionally broad. For momenta close 
to the FS, spectral weight is shifted away from 6 = 0, producing a weak secondary maximum, 
which might suggest a "shadow band" [Q, although the real part of the denominator of the 
Green's function still has only a single zero. These features produce a weak pseudogap in the 
total DOS, even though the FS spectral functions are single-peaked 0. For momenta on or 
near the FS, the real part of the self-energy has an anomalous positive slope near e = and 
the quasiparticle lifetime has a local minimum there. A simple analytic model shows that 
these anomalies result from the coupling of quasiparticles to AFM spin fluctuations, without 
a phase transition to AFM order. The process is anisotropic, beginning at the X-point and 
spreading over the FS. 

The FEA for the self-energy in a paramagnetic state of the Hubbard model is 



S(r, r) = [xphir, r) + T,,(r, r) + T..(r, r)] ^(r, r) 

+U%p{r,T)G{-r,-T), (1) 

where Xp/i(i"7 ''") = ~G{r, r) G{—r, — r) and Tg-o-, Tpp, and Tpp are propagators for spin fluctu- 
ations, density fluctauations, and singlet pair fluctuations. For example, the spin-fluctuation 
propagator is 

T^^{q,u}^)= - -. (2) 

2 1-Uxph{(l,uj^) 

The Green's function is obtained from Dyson's equation, 

G(k, Bn)'' = G,(k, Bny' - S(k, (3) 



where Go(k, is the non-interacting Green's function, and the self-consistent S is found 



by iteration. We have performed calculations on 128 x 128 lattices of k-points with 
typically 512 Matsubara frequencies using a new massively parallel algorithm that treats 
high frequency contributions exactly . This greatly improved treatment of high-frequency 



information allows us to obtain what we believe are the first meaningful results for the FEA 
spectral functions of the Hubbard model at half-filling. Numerical analytic continuation of 
S to the real- frequency axis is accomplished using Fade approximants [1T3| . 

The extinction of sharp quasiparticles at the X-point is evident from the spectral func- 
tions in Fig. |I]. For U = 1.0 and T = 0.1, the FEA solution shows a sharp quasiparticle at 
e = 0, the Fermi energy (henceforth all energies are in units of t). While an increase in the 
interaction strength and thus an increase in the quasiparticle-quasiparticle scattering rate 
leads generally to a reduction in spectral weight at e = and a broadening of the spectral 
function, the effect shown for U = 2.3 in Fig. |l| is much larger than expected. 

To see this, focus on the self-energies at the X-point shown in Fig. 0. Scaling of the 
U = 1 result by U"^ leads to an Im S comparable to the U = 2.3 result at high frequency, 
but smaller than the true self-energy for e by a factor of ~ 2.5. The FEA self-energy 
at low energy is inconsistent with Fermi liquid theory: an anomalous maximum is clearly 
evident at e = in |ImS|, and Re S at e = has a positive slope which is inconsistent with 
the interpretation of (1 — dReT,{kp, e) / de |e=o)~^ as the quasiparticle pole weight. For the 
same U a self-consistent calculation with only the second-order skeleton diagram yields a 
sharp spectral function and self-energies without these anomalies. 

The dramatic loss of spectral weight at the Fermi surface is reflected in the formation of 
a weak pseudogap in the total density of states as a function of temperature as shown for 
U = 1.57 in Fig. ^ The anomalies in S first appear at the van Hove critical points, and 
spread as T decreases; for a range of U and T, the anomalous excitations near the X-points 
coexist with more conventional quasiparticles on the rest of the Fermi surface. 

The redistribution of spectral weight away from the Fermi energy is a consequence of 
strong antiferromagnetic spin correlations, signaled by the growth of sharp structure in 



the spin-fluctuation T-matrix To-o- when UxphiQ,0) approaches unity at Q = (vr, vr). In a 
Hartree-Fock theory, UxphiQ,^) reaches unity at a sufficiently low temperature, signaling 
a transition to AFM order. Within the FEA, the self-consistently determined quasiparticle 
lifetime regulates the growth of f/Xp/i(Q,0). As shown in Fig. ^, U Xph{Q,0) closely ap- 
proaches unity with decreasing temperature at fixed U, and with increasing U for fixed T. 
As the temperature is decreased for fixed U a sharply peaked structure evolves in To-o-(q, ujm) 
within |Ag| ^ O.OSvr of Q (for U = 2.7 and T = 0.10) and, to an excellent approximation, 
at Um = 0. For real frequencies, Tfj„{Q,Lj) has sharply peaked structure as shown in Fig. 
^. At higher temperatures, as shown for T = 0.15, temperature and correlations broaden 
the peak in UxphiQ, 0) and hence in Re Tg-o- (see Fig. | inset) sufficiently that a region of 
positive slope in Re S(kp', e) at e = is not observed for any U and only a weak maximum 
at ?7 ~ 4.8 is evident in UxphiQ, 0) as a function of U. 

At T = 0.1, further increases in U lead to neither a positive slope in Re E nor sharp 
structure in Im S at £ ^ 0, perhaps because the large quasiparticle width 'smears' the Fermi 
surface and reduces the effect of nesting and van Hove critical points. However, UxphiQ,, 0) 
continues to approach unity monotonically up to the largest U (~ 3.6) for which we have 
solutions at T = 0.1. Hence it is the sharpness of UxphiQ,^) (and Tg-o-) and not simply its 
size that leads to the observed anomalies. 

The entropy S and AFM susceptibility Xafm [0] as a function of temperature are 
shown in Fig. ^ for U = 1.57. At high temperatures, S{T) tracks the U = entropy. As 
the pseudogap opens the entropy turns down, signaling the loss of quasiparticle states. For 
T ^ 0.05, Xafm shows only a modest enhancement and the paramagnetic state remains 
stable, as required in 2D. The nearly singular behavior of To-o-(Q,0) is offset by the vertex 
corrections included in a fully conserving description, as suggested by previous calculations 
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To understand better the self-consistent self-energy, we appeal to a simple model moti- 
vated by a calculation of the ffuctuation conductivity above Tc in superconductors |TB[. The 



anomalies in E are always accompanied by a To-o-(q, uJm) that is strongly peaked near q = Q 



and Um = [|T^. Assuming that the spin-fluctuation contribution to S dominates, 

E(k,£„) ^G(k+Q,e„)t,p, (4) 

where tsp is proportional to the weight of the peak in T„„ within a reciprocal correlation 
length ^"^ of Q, 

tsp = — — T^^iQ + q, u;^ = 0). (5) 

|q|<i/« 

Eqs. (H) and (|^), together with Dyson's equation, can be solved directly. For e ^ the slope 
of Re ^(e) is 1/2. In the full calculation, the slope is also positive when To-o- is sharply peaked, 
but generally differs from 1/2 due to other contributions to the self-energy not included in 
Eq. ^. We observe that for slopes greater than unity, multiple quasiparticle peaks can appear 
in the spectral function, corresponding to multiple solutions of e — — Re S(k, £:) = for 
a given wavevector. For T > 0.05 we have not observed such a splitting of the band, 
in contrast to at least two non-self-consistent calculations describing the effects of strong 
antiferromagnetic correlations on quasiparticle properties |]18| , |19| . 



The single-particle spectral function for the model is nonzero only for < < Msp, 
and given by 
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A{k,e) = Jtsp-e.e+lA, (6) 



where £_ = £ — Ck and = e — Ck+q = e + Ck. This spectral function shares two important 
features with the full calculation. First, spectral weight at ek is shifted away from e = 0, 
with the non-interacting delta function peak becoming a square root singularity truncated 
at ek. This transfer of spectral weight from e = leads to the formation of a pseudogap 
in the total density of states similar to the FEA pseudogap shown in Fig. ^ Second, an 
additional peak forms at ek+q as expected from dynamical coupling to antiferromagnetic 
spin fluctuations. We emphasize that ek+q — Ck — Re S(k, ek+q) 7^ so that this second peak 
is not a second quasiparticle solution. Spectral functions for isp = 0.05, a value consistent 
with To-o-(q, tUm) for U = 2.7 and T = 0.10, are shown together with those from the full 
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FEA calculation in Fig. |T]. Given the simplicity of the approximation, the agreement is 
good. Results obtained from analytic continuation of quantum Monte Carlo data at low 
temperature (T = 0.1) and moderate coupling {U = 4) appear consistent with these general 
features |^ . 

In summary, the FEA yields thermally broadened quasiparticles at high T for any U ^ 
8t. With decreasing temperature and for a modest U, a sharply-peaked spin fluctuation 
propagator develops, whose coupling to the evolving quasiparticle excitations leads to a 
radical breakdown of Fermi liquid theory and to a weak pseudogap in the single-particle 
DOS, in the absence of AFM order or of a gap in the spectral functions on the Fermi 
surface. 
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FIGURES 

FIG. 1. (left) Spectral functions calculated in the FEA at T = 0.10, showing a dramatic 

reduction of spectral weight at the X-point (on the Fermi surface) with a modest increase in U 
between T = 0.10 and U = 1.0 (dotted) and 2.3 (dashed). Also shown for U = 2.3 (sohd) is 
a calculation with only the second order skeleton diagram, (right) FEA spectral functions for 
k = (0.89l7r, 0) (solid) and k = (vr, 0) (dashed) shown in comparison with those for the simple spin 
fluctuation model (long-dash-short-dash and dotted respectively) for tgp = 0.05. 

FIG. 2. The real and imaginary parts of the FEA self-energy at the X-point for U = 1.0 
(dashed) and U = 2.3 (solid). Note the 'inverted peak' at £ = in the imaginary part and the 
positive slope at low energy in the real part. Also shown is the self-energy from a self-consistent 
calculation with only the second order skeleton diagram for U = 2.3 (dash-dotted) which is quali- 
tatively similar to the FEA for = 1.0. 

FIG. 3. (top) The density of states as a function of energy for U = 1.57 and T = 0.16 (dashed), 
0.14 (dash-dot), 0.07 (dotted) and 0.05 (sohd) showing the evolution of a weak pscudogap with 
decreasing temperature, (bottom) The self-consistent AFM spin susceptibility (diamonds) and 
entropy (•) as a function of temperature for U = 1.57. Also shown is the entropy for for = 
(dashed). Note that for clarity, xafm{T) has been plotted as XAFMiT)/4xAFM{0-05)- There is 
no sign of an SDW with q = (tt, tt) in the full spin-susceptibility. 

FIG. 4. The self-consistent Ux{Q,0) (top) for U = 1.57 as a function of T and (bottom) for 
T = 0.05 (squares), 0.10 (diamonds), and 0.15 (•) as a function of U. The inset shows analytic 
continuations of the self-consistent T-matrix Re Tf^a{Q,oj) for T = 0.10 and U = 2.7 (solid) and 
for r = 0.15 and U = 4.8 (xlO and dashed). The former shows a pseudogap in the density of 
states but the latter does not. 
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